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ABSTRACT 

We present a novel scenario for globular cluster (GC) formation, where the ultra- 
violet (UV) background radiation effectively works so as to produce compact star 
clusters. Recent observations on the age distributions of GCs indicate that many GCs 
formed even after the cosmic reionization epoch. This implies that a significant frac- 
tion of GCs formed in UV background radiation fields. Also, the star formation in an 
early-generation of subgalactic objects may be affected by strong UV radiation from 
preformed massive stars, e.g., Population III stars. Here, we explore the formation of 
GCs in UV radiation fields. For this purpose, we calculate baryon and dark matter 
(DM) dynamics in spherical symmetry, incorporating the self-shielding effects by solv- 
ing the radiative transfer of UV radiation. In addition, we prescribe the star formation 
in cooled gas components and pursue the dynamics of formed stars. As a result, we 
find that the evolution of subgalactic objects in UV background radiation are sepa- 
rated into three types, that is, (1) prompt star formation, where less massive clouds 
(~ I0 5 ~ 8 Mq) are promptly self-shielded and undergo star formation, (2) delayed star 
formation, where photoionized massive clouds (> 10 8 M Q ) collapse despite high ther- 
mal pressure and are eventually self-shielded to form stars in a delayed fashion, and 
(3) supersonic infall, where photoionized less massive clouds (~ 10 5_8 Mq) contract 
with supersonic infall velocity and are self-shielded when a compact core forms. In 
particular, the type (3) is a novel type found in the present simulations, and eventu- 
ally produces a very compact star cluster. The resultant mass-to-light ratios, half-mass 
radii, and velocity dispersions for the three types are compared to the observations of 
GCs, dwarf spheroidals (dSphs), and ultra-compact dwarfs (UCDs). It turns out that 
the properties of star clusters resulting from supersonic infall match well with those of 
observed GCs, whereas the other two types are distinct from GCs. Hence, we conclude 
that supersonic infall in a UV background is a promising mechanism to form GCs. 

Key words: galaxies: formation - globular clusters: general - galaxies: dwarf - radia- 
tive transfer - hydrodynamics 



1 INTRODUCTION 

Globular clusters (GCs) are characterized by their compact- 
ness, old age, and low metallicity. Although the mass range 
of GCs (1O 5_6 M0) is similar to that of dwarf spheroidal 
galaxies (dSphs), GCs have much higher stellar velocity dis- 
persions than those of dSphs. In addition, the velocity dis- 
persions in GCs correlate with their luminosities in a def- 
initely different fashion from that in dSphs. Furthermore, 
the survey of Fornax galaxy cluster has revealed a new 
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class of dwarf galaxies, that is, ultra compact dwarf galax- 
ies (UCDs), which are more compact t han normal dwarf 
galax ies but clearly distinct from GCs (|Drinkwater et at] 
l2003| h According to the virial theorem, the velocity dis- 
persion is expressed as a* oc -J GM/rh, where M and 
rh are the total mass and the half-mass radius of a sys- 
tem. The velocity dispersion is then related to the lumi- 
nosity by o- t oc L 1 / 3 , assuming M oc r\ and L oc M. 
Obser ved GCs, however, show the correlation like cr* oc 
L 1/2 (|Diorgovski et al.|[l997l : iMcLaughlin fc van der Marell 
2005). This correlation can be derived if GC radii are al- 
most independent of mass (| Ashman fc ZepJ^OOll ). Interest- 
ingly, such a correlation is found not only in Galactic GC s 
but also in GCs in external galaxies fHas egan et al.ll2005l ). 
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Also, the a* — L re lation of UCDs obvious ly deviates from 
the virial theorem (|Drinkwater et all [2003). As for dSphs, 
some follow cr* oc L 1//3 , but the majority of dSphs show 
higher velocity dispersions than expec t ed with the virial the- 
orem jMatedll99a: iMatin et al.ll200rj; ISimon fc Gehall2007l ; 



iMcConnachie fc Irwinl 120061 : iMaiewski et all 120071 ). These 
imply that GCs as well as UCDs and many dSphs did not 
form through simple gravitational collapse and virialization. 

Various scenarios for the formation GCs have 
been hitherto cons i dered by many authors (e.g., 
Peebles fc Dickd 1 19681: iFall fc Reesl 1 19851: [Ashman fc Zepj 
19921: iMurrav fc Linl 1 19921: iBromm fc Clarke! 120021 : 
Kravtsov fc GnedinJ 120051 : ISaitoh et alj 120061 ). They 
are mainly divided into three categories, based on the 
formation epoch of GCs relative to that of a host galaxy. 
The first one is a merger scenario, in which GCs form in the 
gas-rich mergers of disk galaxies l| Ashman fc Zepjl 19921 ). In 
this scenario, interstellar gas is highly compressed by shocks 
induced by the mergers, and very compact star clusters are 
expected to form. Such events likely occur at a later stage 
of galaxy formation, therefore formed star clusters would 
be young. This scenario can explain young star clusters 
observed in merging or starbursting galaxies. However, it 
might be difficult to explain old GC populations. Moreover, 
the metallicity of star clusters is expected to be the same 
as that of merging galaxies, which may be higher than the 
GC metallicity of Z ~ 10" 4 - 10" 2 Z Q . 

In the second scenario, GCs form a t the e arly stage of 
the host galaxy evolution. IFall fc Reed ljl985l ) firstly have 
shown the possibility of GCs formation caused by thermal 
instability, in which proto-globular clusters are formed from 
dense gas clouds with T ~ 10 4 K confined by diffuse hot gas 
with T ~ 10 6 K. In this case, the Jeans mass for such dense 
clouds roughly corresponds to the typical mass scale for ob- 
served GCs (1O 5_6 M0). In this model, the clouds should be 
metal-poor to maintain the cloud tem perature as T ~ 10 4 K . 
A similar idea have been explored bv lMurrav fc Linl (|l992l ). 
in which the temperature of dense clouds is maintained 
as T ~ 10 4 K by photoheating. iKravtsov fc GnedinJ (|2005l ) 
have carried out high resolution simulations, and shown that 
GCs can form in high-density cores of giant molecular clouds 
in dense gaseous disk of galaxies. They also have found that 
the resultant ratio of the total GC mass to the baryonic 
mass of the host galaxy is as small as » 10~ 4 , as indicated 
by observations (|McLaughlinl 11999). 

The last one is the pregalactic formation scenario that 
has been widely explored, in which GCs form in dark mat- 
ter (DM) minihaloes before they infall into galactic haloes. 
This is a plausible scenario for the formation of low-mass 
star clusters, and it can account for many properties of GCs 
such as low m etallicity and old age comparable to the cos- 
mic a ge (e.g., Bromm & Clarkd I2OO2I ; iMashchenko fc Silld 
120051 : ISaitoh et al.ll2006l ). However. proto-GCs in this sce- 
nario are expected to be dominated by dark matter (DM), 
and possess the mass-to-light ratios of M/L ~ 10. There- 
fore, it does not provide a direct explanation for observed 
GCs that have mass-t o-light ratios as low as M/L ~ 1 
l|Prvor fc Mevlanlll993l) Recen t simulations of galaxy for- 
mation by Saitoh et alj (|2006l ) have indicated that diffuse 
DM haloes are selectively stripped by tidal interactions with 
their host galaxy, and baryon-dominated star clusters are 
naturally explained. However, the spatial resolution of the 



simulations cannot allow us to discuss the internal structure 
of proto-globular clouds. Therefore, the physical reason for 
condensations of the proto-globular clouds is still uncertain. 
The stripping of DM haloes also has been s tudied using N- 
body calculations by IMashchenko fc Silld (|2005l ). where a 
hybrid cluster composed of stars and DM is simulated in 
the external tidal field of the host galaxy. They assumed 
initially a compact baryon core embedded in a diffuse DM 
halo, and explored the dynamical evolution using iV-body 
simulations. They have found that the surface brightness 
profiles and mass-to-light ratios for the evolved hybrid GCs 
are consistent with those of observed GCs. However, it has 
not been clarified yet how such compact stellar systems form 
at first. Since the virial radii of collapsed objects at z ~ 10 
are several hundred parsecs in the mass range of 1O 5_6 M0, 
that are ~ 10 — 100 times larger than the sizes of observed 
GCs. Hence it is not easy to form baryon-dominated star 
clusters. An additional physical mechanism seems to be req- 
uisite to produce a compact stellar core. Moreover, for all 
the above scenarios, the physical mechanism to explain the 
cr* — L relation of GCs has not been elucidated. 

The ages of Galactic GCs show a broad dis- 
tribution with the m e an age of about 12.3Gyr 
|Puzia. Perrett fc Bridge! 120051 ). A significant fraction 
of GCs seem to form after the cosmic reionization epoch 
z T ~ 11, which is inferred by Wilkinson Microwave 
Anis otropic Probe (WM AP) three year data (|Page et al.l 
120071 ) and five year data (|Komatsu et al.ll2009l ). Therefore, 
not a few GCs are thought to form in UV background 
radiation fields. Also, the formation of GCs itself might 
contribute much to the reion ization of the universe at 
redshift z ~ 6 l|Ricottill200l ). Furthermore, even before 
the reionization, the star formation in an early-generation 
of galaxies can be affected by strong UV radiation from 
preformed massive stars including Population III (Pop 
III) stars. Thus, it is of great significance to consider 
the formation of GCs in UV radiation fields. When a 
gas cloud is photo-ionized, the cloud is heated up to 
~ 1 4 K (U mcmura fc Ikeuchil 11984 iThoul fc Weinberg! 
1996). In addition, soft UV radiation (Lyman- Werner 
band) destructs H2 molecules, which is the main coolant 
below 10 4 K in low metallicity g as (Steelier & Wi lliams! 



below 10 K m low metallicity g as (btcchcr & William; 
Il967l : lHaiman. Rees. fc Loebl 119971 : IKitavama et all l200ll ) 



Then, the formation of low mass objects with the virial 
temperature less than 10 4 K is suppressed. Moreover, no 
star is expected to form in ionized gas. Hence, to form 
star clusters, the gas cloud should be self-shielded from 
UV background radiation l|Taiiri fc Umemural Il998l ). To 
treat the self-shiel ding, we shoul d solve the radiative 
transfer properly. IKitavama et alj l|200ll ) explored the 
formation of dwarf galaxies within UV radiation fields, 
using spherical symmetric hydrodynamics coupled with the 
radiative transfer. With the radiation hydrodynamic (RHD) 
simulations, they showed that dwarf galaxies can form even 
after the reionizat i on, if the self-shielding effectively works. 
ISusa fc Umemural l|2004) also studied this issue by three- 
dimensional RHD simulations including formation of stars 
and the stellar dynamics. In their simulations, they showed 
that a dwarf galaxy with the stellar mass of » 10 6 A'/ Q is 
able to form even in UV radiation. The half-mass radius of 
the stellar system is ~ lOOpc, which is corresponding to the 
sizes of dSphs, but greater by about one-order magnitude 
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than the observed GC size. Hence, a further mechanism is 
required to produce much more compact star clusters in 
UV radiation. 

In this paper, we explore a possibility that a superson- 
ically contracting cloud produces a compact star cluster in 
UV radiation. A gas cloud with its infall velocity exceeding 
the sound speed of ionized gas (T ~ 10 4 K) can keep con- 
tracting, even if the gas cloud is fully ionized. The contract- 
ing cloud is eventually shielded from both ionizing photons 
and H2 dissociating photons, and then the cloud can cool 
via H2 cooling, generating a star cluster. The star cluster 
formed through such a process becomes very compact due 
to strong energy dissipation. 

This paper is organized as follows. In §2, the simula- 
tion code and numerical procedure are described. The re- 
sults of simulations are presented in §3, and it is shown that 
the evolution of subgalactic objects are typically separated 
into three branches. The comparison with the observations 
of GCs, dSphs, and UCDs is given in §4. The specific fre- 
quencies of GCs, required UV intensity, tidal effects by host 
galaxy, and internal feedbacks are discussed in §5. §6 is de- 
voted to the conclusions. Throughout this paper, we assume 
a A-CDM cosmology with the matter density £Im = 0.3, 
the cosmological constant Q\ = 0.7, the Hubble constant 
h — 0.7, and the baryon density Qb = 0.05. 



2 SIMULATIONS 
2.1 Numerical Scheme 

We use a radi at ion-h ydrodynamic scheme developed by 
iKitavama et alj (|200ll ). to solve hydrodynamics coupled 
with radiative transfer. This scheme is based on the second- 
order Lagrangian finite-difference method in spherical sym- 
metry. In this scheme, we treat self-consistently gravita- 
tional force of dark matter and baryon, hydrodynamics, non- 
equilibrium chemistry of primodial gas including H2, and the 
radiative transfer of ionizing photons. For self-shielding by 
H2 against photo-dissociating ph otons, we employ the sel f- 
shielding function introduced by iDraine fc Bertoldil (|l996T) . 
[See also IKitavama et alj l|200ll ) for the details of scheme.] 
Initially, the number of shells is Nb ~ 600 for gas compo- 
nent and iVd = 10000 for DM component. Shocks are treated 
with artificial viscosity. In order to investigate whether the 
results are sensitive to the treatment of artificial viscosity, 
we perform the simulations with stronger artificial viscosity 
(10 times as strong as the fiducial one) or weaker artificial 
viscosity (0.1 times as weak as the fiducial one). As a result, 
we find that the changes of results are only a few percent. 
Thus, we conclude that the present results are not sensitive 
to the treatment of artificial viscosity. 

We also consider the formation of stars from cooled 
gas component and pursue the relaxation process of formed 
stars. We take the following conditions as the star formation 
criteria: 



dt 



< 0. 



(4) 



T g < 2000A', 



v r < 
dp. 



dt 



>0, 



(1) 
(2) 

(3) 



Here T g , v T , and p g are the temperature, the infall velocity, 
and the density of gas shell. Ionized gas cannot cool below 
10 4 K, if satisfying the condition (3). Hence, the gas must be 
shielded against UV radiation to satisfy the condition (1). In 
that sense, the condition (1) is a most important criterion. 
However, the threshold temperature 2000K is not so crucial. 
When we changed this threshold temperature from 1000K 
to 5000K, we found no significant difference in the main 
results. We assume that if a gas shell satisfies all of the above 
conditions, the gas shell becomes a stellar shell immediately, 
since the timescale of star formation is basically the free-fall 
time in a high density region, which is considerably shorter 
than the dynamical timescale of the whole cloud. 

The basic equation for the stellar dynamics is given by 

d 2 r s _ GM tot (< r s ) 
dt 2 ~ rf 



(5) 



where r s and iWtot(< r s ) are the radius of stellar shell and 
the total mass (including baryon and DM components) in- 
side r s , respectively. For stellar shells, we allow them to cross 
each other. We also calculate the mass weighted velocity dis- 
persion, which is represented by 



2 



(6) 



where N*, dm*, n, v* xx and A7* ; tot are the number of stellar 
shells, mass of n-th stellar shell, infall velocity of n-th shell 
and the total stellar mass, respectively. 

2.2 Setup 

The initial density distribution of a cloud consi s ting o f gas 
and DM is the same as that in IKitavama et al.l (|200lh . We 
calculate the evolution of the cloud from its linear overden- 
sity stage. We concentrate on low-mass clouds with the ini- 
tial baryonic mass of 10 5 M Q < Af b ,in < 10 8 M Q collapsing 
at high redshifts 3 < z c < 20. Since we define z c as the epoch 
at which a cloud as a whole collapses, the collapse redshift of 
the central high-density regions z c o is earlier than z c . These 
two redshifts are related by 1 + z c o = 2.7(1 + z c ). When 
z ^ zuv, the cloud is exposed to external UV radiation. 
Here, we assume the constant UV intensity specified by 72 1 , 
where 72i is the intensity at the Lyman limit frequency of 
hydrogen i>l in units of 10 -21 erg cm _2 s _1 Hz _1 str _1 . 

The cosmic r e ionization epoch is inf erred to be z r ~ 11 
l|Page et alj 120071 ; iKomatsu et all 120091 ). The reionization 
of intergalactic matter can start at earlier epochs in an 
inhomogeneous fashio n (Na kamoto. Umemura fc Susal 
120011 ; ICiardi et alj 12001). Although the ionization sources 
are still uncertain, the reionization by Pop III stars 

has bee n co nsidered as one of p lausible possibilities 

ilGnedinll200"ol ; ICiardi et al.ll200ll ; ICenll2003l ; ISokasian et ail 
2004). Here, we assume Pop Ill-type ionization sources. 
Previous studies have sho wn that Pop III stars are as 



massive as 



lOOAfp (e .g. , lAbel. Bryan fc Normanl |2000|; 



Bromm. Coppi fc Larsor] |2002| ; iNakamura fc Umemural 
200 it lYoshida et alj|2006l ). Then, the blackbody radiation 
with effective temperature of T e g ~ 10 5 K is emitted from a 
Pop III star. The numbe r of ionizing ph otons emitted per 
second is N ~ 10 50 s _1 (jSchaereij [2003 1. Since a Pop III 
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object is as massive as 10 6 M Q (e.g., IOmukai fc Nishi|[l99i 
lYoshida et al] 120031 '), almost all ionizing photons emitted 
from a Pop III star c an escape from the Pop III object 
|Kitavama et alj |2004 ) . Qualitatively, the virial radius of 
the Pop III object is given by 



1.60 x 10' 



( Mi 



\10 6 M Q 



1/3 



(—) 

\l + zj 



pc, 



(7) 



where Mf is the total mass of the Pop III object and z c is its 
collapse redishift. On the other hand, the Stromgren sphere, 
which is estimated for the cosmic mean density, is extended 
to the radius 



R s = 5.45 x 10 3 



N 

1Q 50 



1/3 



20 
1 + J 



pc. 



(8) 



If the radiation flux from a Pop III star is translated into an 
averaged intensity /„, we have 

2 

J„ = fl„(Teff) ( — ) . (9) 



,(T eff ) (f) 



Under the assumption that one Pop III star is born in a halo, 
the UV intensity is evaluated to be from I21 ~ 10 3 at the 
virial radius of PopIII halo to I21 ~ 10~ 3 at the Stromgren 
radius. Hence, we investigate the range of 10~ 3 ^ I21 ^ 10 3 . 
As for UV irradiation epoch, we investigate three cases as 
zuv=15, 20, or 25. 



Prompt star formation 
(strong evaporation) 




Star Formation 
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Delayed star formation 
(weak evaporation) 




Star Formation Relaxation 
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3 RESULTS OF SIMULATIONS 

3.1 Hydrodynamics and star formation 

If a cloud is not irradiated by UV radiation or self-shielded 
from a UV background, it cools down to several 100K by 
H2 cooling. Then, a cloud more massive than the Jeans 
mass of 100K, Mj(100K) ~ 10 5 Af Q , can collapse to form 
stars. When a cloud is photoionized, the gas temperature is 
raised up to ~ 10 4 K, and then the Jeans mass increases to 
Mj(10 4 K) ~ 10 s M e . If the cloud mass is below A/j(10 4 K), 
the ionized gas cannot be confined to a DM halo, but evap- 
orated. Nevertheless, if the infall velocity exceeds the sound 
speed of 10 4 K, the ionized gas can collapse even when the 
total mass is below Mj (10 4 K). Thus, the hydrodynamic evo- 
lution of gas cloud and subsequent star formation are sen- 
sitively dependent on the total mass of cloud, the strength 
of self-shielding, and the infall velocity. In the following, we 
show the typical numerical results, which are basically cat- 
egorized into three types of evolutionary branches. 

3.1.1 Prompt star formation (strong evaporation) 

This is the case of the cloud mass between Mj (100K) 
and Mj (10 4 K). When external UV radiation irradiates the 
cloud, the inner regions at r < i? s hi c id are promptly self- 
shielded to form stars. But, the outer envelope of cloud is 
photoionized and evaporated due to the enhanced thermal 
pressure, because M < Afj(10 4 K). The schematic view is 
shown in Fig. [1] (top panel). The cloud evolution is shown in 
Fig. [2] where the mass shell evolution is shown around the 
self-shielding radius i? s hi c id- After the cloud is irradiated by 
external UV, a shell inside i? s hi c id can collapse, because the 
self-shielded mass is higher than Mj(100K). Eventually, the 



Supersonic infall 
(mild evaporation) 




M,(T- 100K) < M m < M,(T- W'K) 
1 ► 



Figure 1. Schematic views for the prompt star formation case 
(top panel), the delayed star formation case (middle panel), and 
the supersonic infall case (bottom panel). In each panel, ar- 
rows indicate infall velocity vectors. Blue circles indicate the self- 
shielding radii iJ s hieldi inside which the gas is impervious to ex- 
ternal UV radiation and therefore can cool by H2 cooling to form 
stars. Red circles indicate the sonic point -R son ici where the infall 
velocity corresponds to the sound speed of 10 4 K. 



shell cools down, forming stars promptly. In this case, the 
star formation criteria is satisfied from inside out to -Rsnioid- 
The duration of star formation is roughly 10-100Myr. The 
duration is longer for higher mass clouds, since a larger part 
of massive clouds can be shielded. On the other hand, the 
regions outside f? s hi e id are ionized and photoheated. Then, 
shells there evaporate shortly. 

As a result, the mass-to-light ratios are expected to 
be higher owing to the mass loss by the photoevapora- 
tion. Such behaviors can be basically understood as the 
formation mechanism of dwarf galaxies in a UV back- 
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Figure 2. Cloud evolution in the case of prompt star formation. 
The physical quantities around the self-shielding radius are shown 
as a function of redshift, with (a) radii, (b) infall velocity, (c) 
hydrogen number density, (d) electron fraction, (e) H2 fraction, 
and (f) gas temperature of the cloud. The parameters here are 
ztjv = 20, Aft, ; n = 10 6 Mq, z c = 7, and I21 = 1. In each panel, 
a thick solid line shows the shell that finally collapses, while a 
thick dashed line shows the shell that finally evaporates. A thin 
vertical solid line indicates the epoch of UV irradiation, ztjv. 
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Figure 3. Cloud evolution in the case of delayed star formation. 
The physical quantities same as those in Fig. [2] are shown. The 
parameters here are ztjv = 20, ; n = 10 s Mq, z c = 7, and 
I21 = 100. The time evolution of the outermost shell that finally 
collapses is shown. 



ground tha t is previously studied in detail by several au- 
thors (e.g., Kitavama et al.l l200ll ; ISusa fc Umemural |2004 
iRicotti. Gnedin fc Shullll200S '~ 



3.1.2 Delayed star formation (weak evaporation) 

This is the case of cloud mass slightly higher than 
Mj (10 4 K). When UV intensity is very weak, the almost all 
regions of cloud are self-shielded and collapse in a similar 
way to the UV-fxee case. However, if external UV is rel- 
atively strong, the bulk of cloud is photoionized. No star 
forms in the ionized regions, since no H2 cooling works there. 




24 22 20 18 
redshift 
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Figure 4. Cloud evolution in the case of supersonic infall. Same 
as Fig.[2] except for M b in = 10 7 M Q and I21 = 100. In panel (b), 
the sound speed of 10 4 K is shown by the horizontal dotted line. 
The time evolution of shells around the sonic point is shown. 



Nonetheless, most of ionized regions do not evaporate but 
collapse, because the cloud is more massive than Mj (10 4 K) 
and therefore the self-gravity overwhelms the pressure of the 
ionized gas. As a result, in the course of cloud contraction, 
the ionized regions are self-shielded due to the increase of 
density. Then, stars form there in a delayed fashion. In this 
case, the star formation can continue for > lOOMyr because 
of the delayed self-shielding. The schematic view is shown 
in Fig. [T] (middle panel). The cloud evolution in a strong 
UV background is shown in Fig. [3] where the outermost 
shell that finally collapses is shown. The shell is ionized af- 
ter ztjv = 20, but it collapses and is self-shielded around 
z = 14, forming a cooled shell. As a result of such de- 
layed self-shielding, a formed star cluster can become more 
compact than a star cluster formed in no UV background, 
and therefore the stellar velocity dispersion of cluster is in- 
creased. It is noted that the effect of delayed star formation 
becomes conspicuous for the mass slightly above Mj (10 4 K), 
since a much higher mass cloud is impervious to photoion- 
ization except for surface thin layer. 



3.1.3 Supersonic infall (mild evaporation) 

This is a novel branch that is found in the present simula- 
tions. The cloud mass is between Mj(100K) and Mj(10 4 K), 
similar to the case of prompt star formation, but the infall 
velocity of the cloud can exceed the sound speed of 10 4 K 
owing to radiative cooling as well as a potential of DM halo. 
When the cloud is exposed to strong UV radiation, ionizing 
photons permeate the cloud deeply, and a large part of cloud 
is photoionized. However, shells collapsing with supersonic 
infall velocities cannot be stopped by photoheating. There- 
fore, they keep contracting, regardless of photoionization. 
Such shells are eventually self-shielded, when the density in- 
creases sufficiently due to the contraction, and then forms a 
compact star cluster. The schematic view of the supersonic 
infall case is shown in Fig. Q] (bottom panel). The evolution 
of shells around the sonic point is shown in Fig. [4] A dot- 
ted line shows an evaporating shell, while a solid line shows 
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Figure 5. Dependence of cloud evolution on the cloud mass, col- 
lapse epoch, UV intensity, and UV irradiation epoch. The collapse 
redshift z c and the type of evolution are shown against the initial 
baryonic mass of cloud, ; n . Panels (a), (b) and (c) show the re- 
sults with ztjv = 20, but for different UV intensity as I21 = 0.01, 
I21 = 1, and I21 = 100, respectively. Panel (d) is the result of 
ztjv = 15 and I21 = 1. Triangles, squares, and circles indicate 
the cases of prompt star formation, delayed star formation, and 
supersonic infall, respectively. Pluses represent complete evapo- 
ration, and asterisks correspond to evolution in no UV. The per- 
centages attached to the symbols represent the fraction of the 
final stellar mass Af* to the initial baryonic mass, M t /M^, in . In 
each panel, a long-dashed and a short-dashed line correspond to 
la and 2a CDM fluctuations. A solid line represents the Jeans 
mass of ionized gas Mj(10 4 K). 



a collapsing shell with supersonic infall velocity. It is clear 
that the supersonic infall shell is photoionized once (panel 
b) , and eventually self-shielded from both ionizing photons 
(panel d) and HVdissociating photons (panel e), and cools 
below 10 4 K (panel f) by H2 cooling. Hence, the gas inside 
the sonic point i? SO nic results in a compact star cluster. 

As a result, the strong contraction before star formation 
enhances the stellar velocity dispersion to a large degree. 
As a matter of importance, in this supersonic infall case, 
the duration of star formation becomes less than lOMyr, in 
contrast to the other two cases. This means almost coeval 
star formation, which is consistent with the fact that GCs 
commonly have a single stellar population. The supersonic 
infall branch appears if the sonic point is larger than the 
self-shielding radius. 



3.2 Parameter Dependence 

The dependence of cloud evolution on the mass, collapse 
epoch, UV intensity, and UV irradiation epoch are summa- 
rized in Fig. E] In each panel, the resultant types of evolution 
are shown in the diagram of the collapse redshift against the 
initial baryonic mass of cloud. Also, la and 2a CDM density 
fluctuation spectra are shown. The UV intensity is changed 
with J21 = 10~ 2 ,1 or 10 2 . In the panels (a), (b) and (c), 
the UV irradiation epoch is assumed to be ztjv = 20, while 
ztjv = 15 in the panel (d). In the simulations with ztjv = 15, 
the objects collapsed before z = 15 correspond to the UV- 
free case. Interestingly, the boundary between the complete 
evaporation and star cluster formation roughly corresponds 
to la CDM fluctuation. 

In each panel, a thin solid line indicates the Jeans mass 
of ionized gas Mj(10 4 K), which is given by 



Above this line, the cloud evolution results in the delayed 
star formation as seen in Fig. [5] Below the Jeans mass 
Mj(10 4 K), clouds with more massive and higher collapse 
epochs lead to the star cluster formation by the super- 
sonic infall . This tendency can be understood by follow- 
ing two reasons: (i) The self-shielding radius is weakly de- 
pendent on the cloud m ass (-R shic id oc M 2 ^ 3 ) shown by 
iTaiiri fc Umemural (| 1998f ) , while the radius where the infall 
velocity exceeding the sound speed is roughly proportional 
to the cloud mass. Hence, a higher mass cloud has larger 
ionized, supersonic regions, (ii) The fraction of supersonic 
infall regions becomes larger according as the cloud con- 
tracts. Thus, earlier collapsing clouds tends to form most 
stars through supersonic infall. 

In ztjv = 15 simulations, the trend is the same as the 
simulations with ztjv = 20, as long as the collapse epoch is 
later than ztjv = 15. On the other hand, in the case of weak 
UV with I21 — 10~ 2 (panel a), the supersonic infall case 
does not appear, because the self-shielding radius is always 
larger than the sonic point. 



3.3 Stellar Dynamics 

Here, we pursue the subsequent stellar dynamics of star 
clusters formed through three mechanisms shown above. In 
Fig. [II] the time sequence of mass distributions is shown for 
prompt star formation, delayed star formation, and super- 
sonic infall cases. Obviously, the mass distribution of a star 
cluster formed through the supersonic infall is different from 
the other two cases. In the supersonic infall case, the stel- 
lar component is predominant in the inner several lOpc re- 
gions. On the other hand, in the other two cases, the system 
is dominated by the DM component in almost all regions. 
These results are understood as follows: In the supersonic 
infall case, the star formation is delayed by the external UV 
radiation, until the cloud contraction eventually causes the 
self-shielding. Then, the strong energy dissipation occurs in 
a compact region, leading to the formation of a dense star 
cluster. As a result, a compact star-dominant system forms. 
On the other hand, in the prompt star formation case, stars 
are born at earlier dynamical phase of initially self-shielded 
regions. Consequently, the energy dissipation is not so strong 
and therefore a more diffuse star cluster forms. In the de- 
layed star formation case, the cloud is self-shielded after it 
contracts to some degree. Hence, the dark matter contribu- 
tion becomes smaller compared to the prompt star forma- 
tion case. But, it does not grow into a totally star-dominant 
system, since the cloud is self-shielded before it becomes 
baryon-dominated. Thus, we conclude that the supersonic 
infall is only the branch that produces star-dominant com- 
pact clusters. 



4 COMPARISON WITH OBSERVATIONS 

4.1 Mass-to-light ratios and half-mass radii 

As shown in §3, the final mass-to-light ratios depend on 
the formation processes. Here, we derive the mass-to-light 
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Figure 6. Mass distributions for typical three cases: prompt star 
formation case (top panel), delayed star formation case (middle 
panel), and supersonic infall case (bottom panel). In each panel, 
time variations of mass distributions are shown by M(R) in units 
of the solar mass, where M(R) is the cumulative mass in the 
regions of r ^ R. Thick lines show stellar components, while thin 
lines do DM components. Solid, short-dashed, and long-dashed 
lines are respectively corresponding to the distributions at 8.9 X 
10 6 yr, 2.89 X 10 7 yr and 9.73 X 10 7 yr after z uv . 



ratios of simulated objects. In this section, we use sim- 
ulation results with ztjv = 20, I21 = 10 _2 ,1, or 10 2 , 
Mb,in = 1O 5 M - 1O 8 M , and z c = 5 or 7, which cor- 
responds to the central collapse redshift of z c o ~ 15.2 or 
20.6, respectively. For each simulated object, the dynamical 
mass is evaluated by the total mass (Md yn = M* + A/dm) 
inside the half-mass radius of the stellar component. As 
for stellar components, we assume M*/Lv = 2, where M* 
is the total stellar mas s and Ly is the V-band luminosity 
jPrvor fc Mevlanlll993l ). This is a typical value when stars 
produced by initial starbursts evolve for 10 Gyr. In Fig. [7] 
the resultant mass-to-light ratios Md yn /Lv of simulated ob- 
jects are shown as a function of absolute V-band magnitude 
My, and they are compared to the observations of glob- 
ular clusters (GCs), dwarf spheroidals (dSphs), and ultra- 
compact dwarfs (UCDs). 

As seen clearly, the models of prompt star formation, de- 
layed star formation, and supersonic infall are distinctively 
separated in this diagram. In particular, the model of su- 
personic infall matches fairly well with bright GCs. It is 
worth noting that even if we change Al,/Lv by a factor of 
2, it does not alter basic results significantly. The models of 
prompt star formation and delayed star formation seem to 
match dSphs and UCDs, respectively in this diagram. How- 
ever, these two models can be affected by internal feedbacks. 
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Figure 7. Mass-to-light ratios as a function of absolute V-band 
magnitude My. Large filled symbols are simulations, while small 
open symbols are observations. Filled triangles, filled squares, 
and filled circles indicate the objects formed through prompt 
star formation, delayed star formation, and supersonic infall, 
respectively. On the other hand, open triangles, open squares, 
and open circles indicate observed dwarf spheroidals (dSphs), 
ultra-compact dwarfs (UCDs), globular clusters (GCs), respec- 
tively. The obser vational data of GCs a re taken fro m the data of 
Milky Way GCs ijPryor fc Mevlanlll993h. M31 GCs jFischer et al.l 
Il993l; iDubath et alj Il996l; iDiorgovski et alj Il997ft . LMC GCs 
I Dubath fc Grillmairl Il997l; iFischer et al.l Il993h, SMC GCs 



Dubath et al.lll992h . and NGC5218 GCs [[Martini fc Hoi 1200 „ 
The d at a of dSphs are t aken fr om iMatec lll998l), iMatin et al. 



d2006h. ISimon fc Gehal d2007tl . iMcConnachie fc Irwinl l|2006l) . 
and | M aiewsk ret**ai* I 11200711 . The UCD data are taken from 
iDrinkwater et alj l|200a ^ 
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Figure 8. Half-mass radii R^ as a function of absolute V-band 
magnitude My- The symbols are the same sense as Fig. [7] 



The effects by the internal feedback are discussed meticu- 
lously in the next section. Anyhow, if the final M»/Lv for 
stellar components is in the range of 1 < Af»/Lv < 4, three 
models are distinctive in this diagram. The mass loss by in- 
ternal feedbacks in prompt star formation and delayed star 
formation would result in an up-shift in Fig. [7] leading to a 
further deviation from GCs. 

Also, in Fig. [S] the resultant half-mass radii of sim- 
ulated objects are shown as a function of absolute V-band 
magnitude Mv, and compared the observations. Again, the 
model of supersonic infall is well concordant with observed 
bright GCs, whereas the other two models are not. Here, a 
remarkable property for the supersonic infall model is that 
the half-mass radius is almost independent of the mass. This 
is due to the fact that a higher mass cloud has larger ionized, 
supersonic region, and eventually the energy dissipation oc- 
curs more strongly. It is worth noting that half-mass radii 
in prompt star formation case are about tens times larger 
than those in supersonic infall case for the same luminosity 
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Figure 9. Velocity dispersions <r» as a function of absolute V- 
band magnitude My. In the left panel, the dependence on UV 
intensity is shown. Asterisks, reverse triangles, diamonds, and 
pentagons represent the results for I21 = 0, I21 = 0.01, I21 = 1, 
and I21 = 100, respectively. The virial relation expected in the 
UV-free case are also shown by a long-dashed, a solid, and a 
short-dashed line, respectively for the collapse epoch of z c =5, 7, 
and 10. In the right panel, the simulation results are compared 
with the observational data, where the meanings of symbols are 
the same as Fig. [8] The best fit relation for the observed GCs is 
shown by a dotted line. 

range. These results imply that UV background radiation 
works significantly to determine the properties of subgalac- 
tic objects. The final mass-to-light ratios and half-mass radii 
are both responsible for the velocity dispersions of objects. 

4.2 Velocity dispersions 

In Fig. [9l the stellar velocity dispersion cr» of simulated ob- 
jects are shown as a function of absolute V-band magnitude 
My. Here, we also show the results with I21 = (UV-free 
case), assuming M*/Ly = 2. In the left panel, the depen- 
dence on UV background intensity is shown for I21 = 0, 
I21 = 0.01, I21 = 1, and I21 = 100. Also, the relation 
(L oc crj) predicted from the virial theorem in the UV-free 
case is shown for the collapse redshifts of z c =5, 7, and 10. 
The objects without UV effects well follow this prediction. 
In the presence of the UV background radiation, systemat- 
ically larger velocity dispersions are predicted for a given 
V-band magnitude. 

In the right panel, the simulation results are compared 
with observational data. The objects formed through su- 
personic infall are distributed near the best fit relation for 
the observed GCs shown by a dotted line. It should be noted 
that only the supersonic infall case can account for GCs with 
high velocity dispersions > 10 km s" 1 . For the supersonic 
infall case, higher mass clouds result in the stronger energy 
dissipation. Consequently, the relation between cr* and Mv 
becomes steeper than that the UV-free case (cr, oc L 1 ^ 3 ). 
The correlation for the supersonic infall case is roughly given 
by cr , oc L 1 / 2 , which is consistent with the observed GCs 
(e.g.. lMcLaughlin fc van der M arcl 2005). Since velocity dis- 
persions are estimated as cr» oc -J GM/R, this relation im- 
plies that the radii of objects is almost regardless of the 
mass. This is consistent with the resultant half-mass radii 
Rh shown in Fig. [8] 



The velocity dispersions for the delayed star formation 
case are systematically larger than the prediction by the sim- 
ple virial theorem, since the energy dissipation is stronger 
than that in the UV-free case. In this diagram, the velocity 
dispersions for supersonic inf ailing case and prompt star for- 
mation case are degenerated at —8 ^ Mv ^ —10. But, the 
origin is different, since the distributions in mass-to-light ra- 
tios and half-mass radii are well separated as shown above. 
Since the outer regions of cloud evaporates in the prompt 
star formation case, the velocity dispersion is determined 
basically by the DM component. 

In Figs. [7] [5] and [§] the supersonic infall does not re- 
produce faint GCs in the range of —6 > Mv > —9- If we 
subdivide the parameter space near the low-mass boundary 
between the models of supersonic infall and prompt star for- 
mation, compact star clusters like faint GCs possibly appear. 
Another possibility is the mass loss through tidal stripping 
by a host galaxy, which results in low-mass GCs. 



5 DISCUSSION 

5.1 Specific frequencies 

It is known that early-type galaxies ha ve higher spe cific fre- 
quency of GCs than late- type galaxies (Harris 1991). More- 
over, the spe cific frequency depends on the luminosity of 
host galaxy l|Forbesl 120051 ; iBekki et all I2006T ). The specific 
frequency Ss is defined as the GC population normalized to 
Mv.host = 15 as 

Sn = N t X 10°- 4 ( M v,ho at +15) i 

where Nt is the total number of GCs in a host galaxy and 
Mv.host is the V-band magnitude of host galaxy. This ten- 
dency could be qualitatively understood by combining the 
present results with the galaxy formati on theory in UV 
backg r ound radiation, wh i ch is s tudied bv lSusa fe Umemural 
( 20001. ISusa &: Umemural (|2000T ) have shown that early-type 
galaxies form more preferentially at high-cr peak regions of 
CDM fluctuations, since the self-shielding works effectively 
to allow the high efficiency of star formation and therefore 
the galaxy formation proceeds in a dissipationless fashion. 
The present simulations have shown that GCs can also form 
preferentially from high-cr fluctuations (> 2cr) as shown in 
Fig. [5] Therefore, it is expected that the number of GCs 
around an early-type galaxy tends to be larger than that in 
a la te-type galax y. 

Im oorc et al l (|2006h have explored the spatial distribu- 
tion of subhalos in host galaxies, by high-resolution Af-body 
simulations. They have found that the spatial distributions 
of subhaloes originating in high-cr peaks (> 2.5cr) are simi- 
lar to those of old metal-poor GCs in the Milky way galaxy. 
Thus, it is likely that objects by the supersonic infall can 
corr espond to the meta l-poor GCs in host galaxies. 

ISaitoh et al.l (I2006I - ) have carried out hydrodynamic sim- 
ulations on the galaxy formation with high spatial resolu- 
tion. It is found that numerous clumps of globular cluster 
mass scale form in a host galaxy. The specific frequency 
iSn can be roughly estimated to be the Sn ~ 20, assuming 
the mass-to-light ratio of M/Lv = 5. On the other hand, 
the speci fic frequency in spiral galaxies is observed to be 
arris! 119911 ). This implies that not all the small 
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mass clumps do not evolve into GCs. Here, we have shown 
that only the supersonic infall branch can lead to the forma- 
tion of GCs. Hence, a portion of small mass clumps evolve 
into GCs. To make quantitative argument, it is necessary 
to perform high resolution simulations on the galaxy forma- 
tion including the radiative transfer effect of UV background 
radiation. This would be a future challenge. 

5.2 Sites of GC formation 

As shown in Fig. [5] the supersonic infall case does not 
appear in the simulations with I21 ^ 0.01. If there is no 
or very weak UV radiation fields, only low density star 
clusters form (top panel in Fig. |6}. To form compact star 
clusters like GCs, strong UV radiation, roughly I21 > 1 
is required. Massive Population III stars are one possibil- 
ity to produce strong UV intensity in high redshift epochs. 
But, the strong UV radiation is relative ly localized in HII 
region around a Population HI halo ( Susa fc Umemural 
2006; lHasegawa. Umemura fc~ Susa 2009). Another possibil- 
ity is b lack hole accretion that c ould be a more powerfu l 
source dRicotti fc Ostrikeri |2004| ; ISusa fc Kitavamal |2000| ). 
ISusa fc Kitavamal (|2000l ) simply estimated the UV intensity 
around AGN of 10 44 erg s _1 as 



, / lOOkpc 

lx {— 



. . -1 n-21 -1 -2 . -1 

x 10 erg s cm str , 



(12) 



where R s is the distance from the black hole. That is to 
say, /2i ~ 1 is achieved at lOOkpc around an AGN. This 
can allow the GC formation in a fairly wide area around the 
host galaxy. 



5.3 Tidal stripping 

It is expected that the evolution of star clusters can be af- 
fected by the tidal force by its host galaxy (i.e. Mashchenko 
& Sills 2005) . The importance of the tidal interaction is esti- 
mated as follows. The balance between the self-gravity and 
the tidal force is given by 



Gm 2 {< n) 2GM host m(< r t )r t 

m 



(13) 



where r t ,m(< r t ),Mhost and R c are the tidal radius, the 
cluster mass inside r t , the host galaxy mass and the position 
of star cluster from the centre of host galaxy, respectively. 
GCs are likely to form in small proto-galaxies in the early 
universe (e.g., iBromm fc Clarke! 120021 ; iMashchenko fc Sills! 
2005). Hence, we assume that il4 ost = lO 9 A/0 and R c = 
0.3— lkpc. If we use the distributions obtained by the present 
simulations, the tidal radii r t are evaluated as lOpc - 50pc 
for the supersonic infall case (bottom panel in Fig. [6j . This 
value is larger than the half-mass radii of GCs shown in 
Fig. [5] This implies that the estimation of mass-to-light ra- 
tios or velocity dispersions for the supersonic infall case is 
not affected, even if the outer diffuse DM halo is stripped 
away by the tidal force. Unfortunately, it is difficult to urgue 
precisely the evolution of star clusters in the tidal fields by 
one-dimensional simulations. Hence, we need to simulate the 
three-dimensional dynamical evolution by means of A-body 
method. We are going to present such A-body simulations 
in a forthcoming paper. 



5.4 Internal feedbacks 

Our calculations do not include the internal feedback, i.e., 
UV radiation from internal massive stars and supernova 
(SN) explosions, which becomes important in the case of 
successive star formation. As mentioned above, the dura- 
tion of star formation in supersonic infall case is typically 
several Myrs. The binding energy of a star cluster produced 
by supersonic infall is in the order of 10 51 erg, which corre- 
sponds to the typical SN energy. Hence, only one massive 
SN can sweep away the interstellar gas. Therefore, the su- 
personic infall results in a single stellar population, which is 
consistent with observed GCs. It is also important to men- 
tion that such a short duration of star formation may also 
be responsible for the chemical homogeneity observed in in- 
dividual GCs. 

On the other hand, much larger stellar metallicity dis- 
persions found in dwarf galaxies may be due to succes- 
sive star formation. In the cases of prompt star forma- 
tion and delayed star formation, the duration of star for- 
mation is 10 — lOOOMyr, and the binding energy roughly 
ranges from 10 51 erg to 10 54 erg. The importance of inter- 
nal feedbacks is likely to depend on the mass scale. In the 
case of Mb, in < 10 7 Mq, the duration of star formation is 
shorter (several lOMyr) as shown in £13.11 since the star 
formation is quickly quenched by the external UV radia- 
tion. In addition, SN explosions as well as stellar UV ra- 
diatio n would significantly suppress subsequent star fo rma- 
tion (|Kitavama et alj|2004 iKitavama fc Yoshidall2005h . In- 
terestingly, a recently discovered "fai nt" dSph Ursa Major I 
shows a single old stellar population (|Okamoto et al. I 2OO8). 
The low-luminous prompt star formation model might cor- 
respond to such "faint" dSphs. On the other hand, in the 
case Of M b ,in > 10 7 M Q , successive star formation is possible 
even if the cloud is irradiated by the external UV radiation. 
Moreover, repeated SN explosions, which might trigger the 
star formation, are possible owing to the fact that the bind- 
ing energy is larger than ~ 10 51 erg. Hence, in this case, the 
present simulations should not be directly compared with 
observations. The final st ates would be determined by inter- 
nal feedbacks as shown in lDong. Lin fc M urray (2003). They 
have studied the effect of the internal feedbacks by UV and 
SN explosions of formed stars as well as the external UV 
feedback on the formation of dSphs, by spherical symmetric 
simulations. They have shown that the star formation in low- 
mass dSphs is severely suppressed by the external UV radi- 
ation, whereas the star formation in high-mass dSphs is self- 
regulated by the internal feedback processes, and the final 
equilibria are determined solely by the internal feedback. In- 
terestingly, their simulations have succeeded in reproducing 
multiple stellar population, which is observed in some dSphs 
in the Local Group (e.g. , Grebell 19971; Ikuta fc Arimotol 



l2002j ; lGrebel fc ; Gallagherll2004 iTolstov et alj|2004i r 

This argument can be also applied in delayed star for- 
mation case. Although the final state is likely to be deter- 
mined by the internal feedbacks, the external UV works to 
delay the onset of star formation. Owing to such an effect, 
the energy dissipation becomes slightly strong, which leads 
to the enhancement of velocity dispersion at a central region 
(mid dle panel of Fig. O - Observed UCDs dDrinkwater et ail 
2003) might be explained by this process (see Fig. [9j . 

While the important processes of internal feedbacks 
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has been studied by iDong. Lin fc Murray! ([2003) in spher- 
ical symmetry, there a re also three-dimensional effects. 
ISusa fc Umemural (|2004l) performed three-dimensional ra- 
diation hydrodynamic simulations with solving radiative 
transfer to investigate the shielding effect of local den- 
sity peaks within a UV background. They have shown 
the local shielding allows the long duration of star for- 
mation even after the reionization. In addition, inter- 
stellar medium is compressed intricately by shocks pro- 
duced by multiple SN explosions, and therefore the 
star formation would be tr iggered by such explosions 
II Mori. Ferrara fc Madaull^Ool ). Moreover, we mention that 
the dynamics of stars w ould be more complex, owing to 
three-dimensional ef fects l|Aarseth. Lin fc Papaloizoulll988l ; 
iMurrav fc Linl 1 1996) . These three-dimensional effects can 
play important roles for the star formation history and 
the dynamics of stars. Thus, we plan to carry out three- 
dimensional RHD simulations including the external radi- 
ation as well as the internal feedback to explore the star 
formation history comprehensively. 



6 SUMMARY 

We have carried out the radiation hydrodynamic simulations 
to explore the possibility that the formation of GCs is in- 
duced by external UV radiation fields. As a result, we have 
found that the supersonic infall enables a low-mass gas cloud 
to form a compact star cluster in the external UV radiation. 
A gas cloud with its infall velocity exceeding the sound speed 
can keep contracting, even if the gas cloud is fully ionized 
without the self-shielding from external UV radiation. The 
contracting cloud is shielded from a UV background when a 
compact cloud core forms, and cools by tb cooling. Conse- 
quently, a compact star cluster forms in a diffuse DM halo. 
We have also calculated the dynamical evolution of stars. It 
is found that resultant mass-to-light ratio, half-mass radius, 
and velocity dispersion of simulated star clusters match well 
with those of observed GCs. Therefore, the supersonic infall 
in a UV background is a promising mechanism to form GCs. 
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